Is the Schwarzschild Metric a Vacuum Solution of 
the Einstein Equation? 

Horace Crater* 
The University of Tennessee Space Institute 

June 13, 2011 

Abstract 

This paper examines the inhomogeneous Einstein equation for a static 
spherically symmetric metric with a source term corresponding to perfect 
fluid with p — —p. By a careful treatment of the equation near the origin 
we find an analytic solution for the metric, dependent on a small param- 
eter e, which can be made arbitrarily close to the Schwarzschild solution 
as e — y and which in that same limit can be viewed as arising from a 
point-like source structure. 

1 Introduction 



This paper examines solutions of the inhomogeneous Einstein equation for static 
spherically symmetric metrics corresponding to a perfect fluid at rest with 
p = —p. Our aim is to develope an analytic solution for the metric that can be 
made arbitrarily close to the Schwarzschild solution but yet retains a nonvanish- 
ing contribution to the source term, unlike the Schwarzschild solution. At the 
same time we find that the second solution with the same source term is arbi- 
trarily close to the de Sitter solution. The Einstein equations for the metric are 
second order and highly nonlinear. This implies that if one has two independent 
solutions, then their linear combination will not be a solution. Nevertheless, as 
we emphasize in the first section, under those special circumstances for the 
source term in the Einstein equation, the metric can be written in terms of the 
logarithm of a potential-like function which satisfies a very simple linear differ- 
ential equation, with two linear independent solutions. One of the potential-like 
functions is the Newtonian gravitational potential which through the logarithm 
we associate with the Schwarzschild solution while the other is a harmonic os- 
cillator potential associated with the de Sitter solution. The contibution of the 
latter to the source term is a constant while that of the Schwarzschild solution 
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has a vanishing contribution to the source term. In order to accomplish our 
aim, we recast the inhomogeneous Einstein equation as the limiting case of a 
closely related inhomogeneous equation dependent on a small parameter e. A 
careful treatment of this problem near the origin leads to source terms with two 
separate non-vanishing contributions 



2 Solutions of the Inhomogeneous Einstein Equa- 
tion 

2.1 Static Spherically Symmetric Solutions of the Inho- 
mogeneous Einstein Equation for Perfect Fluid with 

V = -p- 

In this section we first review the static spherically symmetric standard solution 
of the vacuum Einstein equation. The equation is second order and has two 
solutions with one being the Schwarzschild solution while the other is a constant. 
We then remind the reader how if one adds a source term corresponding to a 
perfect fluid at rest with p = — p a second nonconstant solution emerges in 
addition to the Schwarzschild one. If the pressure and density are constants 
the second solution is one found originally by dc Sitter for the Einstein equation 
with a cosmological constant [3]- [5] and is a static form of the time dependent 
one used in models of inflation and dark energy in modern cosmology. 

For a spherically symmetric solution one chooses the coordinates coordinates 

x° = t, 
x = r, 



x 2 = 



(1) 



In a vacuum with static conditions as well as spherical symmetry, we use Dirac's 
exponential parametrization of the metric, pQ, 

dr 2 = e 2v ^dt 2 - e 2X ^dr 2 - r 2 (d9 2 + sin 2 Odcf), 

9oo = -e 2 ^ = l/g m 

gu = e 2 ^)=l/ 9 11 , 

522 - r 2 = l/g 22 7 

g 33 = r 2 sm 2 = l/g 33 . 

g^u = 0, p, ^ v. (2) 



With 
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the only nonzero T's are [I] 



r 1 

1 00 




1 10 






1 ii 


= A', 


r 2 

1 12 






1 22 


= -re~ 2 \ 


r 3 
1 23 


= cot 6, 




1 33 


= -rsin 2 6»e" 2A , 


r 2 

1 33 


= — sin 8 cos 6*. 


(4) 



With 

= r„ A , CT -r w ,j +r£ CT r KA — r£ K r CTA , (5) 

the diagonal elements of the Ricci tensor are 

i?oo = Iv + AV-^—Je^, 

r 

R22 = (1 + rv' — rA')e~ 2A — 1, 

i?33 = i?2 2 sin 2 0. (6) 
From this we have that the scalar curvature is 

R = Q^R^u = - (-v" + AV - v 12 - ^] e- 2A + L" - AV + */ 2 - ^] e- 2A 



2(1 + ri/' - rA')e~ 2A - 2 



4A'-4z/ 2 A _ 9i 2 



2*/' - 2AV + 2*/ 2 - 7 W + -1 e- 2A - ^ (7) 



For our model for T MJ/ we take that of a perfect fluid, [2] 

T^z, = pff^ + {p + p)u,j,u v , 
g^vTu" = -1. (8) 

Consider the fluid at rest, for which u = and so 

5oo u 02 = -1, 

«o = ffoou = .goo(-.9oo) _1/2 , 

«o = ^500- (9) 

Thus, the only nonzero elements of are 



Too 


= P9oo 


Tn 


= pgn, 


T22 


= V922, 


T33 


= P933- 



(10) 
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Now, the Einstein equations G^ v = T^ v become 
G o = -Roo - 2- 900 ^ = ~9ooP 

Goo = i-v" + AV - v' 2 - —) e 2 »- 2X +(v"- AV + v n - 2X ' ~ ^ + \ 
\ r ) \ r r A 

9\' 1 \ p 2v 
' ■ ' A e 2 ^ 2A - %■ = e 2 »p, 



pr 2 



_ 2A / 2A' 1 \ 1 

P = e + 

and 

Gn = Ru - ^9uR = pgu, 



s 2A = i/" - AV + */ 2 - — - \e 2 \ ( 2v" 2AV + 2^ - 4A ' ~ ^' + % 

r 2 \ r r z 



p = e~ 2X 



2u' 1 \ 1 



and 

G22 = -R22 — T^2lR = V912-, 



= (1 + rv 1 - r\')e- 2X - 1 - \r 2 ( f 2*/' - 2A V + 2u' 2 - 4A ' 4l/ + ^ 

2 \ r r z J 



p = (- - -)e~ 2X - dv" - AV + v' 2 - 2A ' ~ 2V ' \ e~ 2X ) 



r r 



v" - AV + v' 2 - — — ] e- 2x 



and the fourth equation, the one for G33, gives nothing new beyond that for 
G22- Hence, the above three simultaneous equations become 



e- 2X 



p = e- 2A 





1 


r 2 ) 






+ 1 








X'-v 1 




r 



p = - \v" - AV + v' z - —-- j e- 2A . (14) 

These are three nonlinear inhomogeneous equations for two unknown functions 
of r. 



For empty space (p = p = 0) these equations become those originally solved 
by Schwarzschild, that is. 



/_AV + /-^) e - 2A = 0. (15) 



Combining the first two equations implies that 



A' = -v\ 

A = -v + X Q (t). (16) 

The third equation then leads 

v" + 2v' 2 + — = 0. (17) 
r 

We parametrize the exponential metric function v by introducing a potential-like 
function V, 

v = iln(l + n 

v > = V- 1 



2 1 + V 

-V"— -V' 2 l —x = -V"— 2u 12 . (18) 

2 1 + V 2 (1 + V) 2 1 + V y ' 

e 2v = (l + V), 
IV 

V" + — = 0. (19) 
r 

This second order equation is an equidimensional one and has the general solu- 
tion of 

V = — + hi. (20) 



Thus, using 

Eq. (fTT|) becomes 



Our metric is thus 



r 



e 2v = l + ^ + h, 
r 

-2i/+2A c 

e 2A = e -2,+2A = ^ (21) 

l + ^ + k 2 
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In order for the metric to become Minkowskian at r — > oo we must have 

fc 2 = = A . (22) 

Matching <?oo to —1 — 20 where for large r where <j) is the Newtonian potential 
—MG/r gives 

ki = -2MG=-r s , (23) 

the Schwarzschild radius. 

Another set of exact solutions is found by assuming that 

P=-P^0. (24) 

Just as with the Schwarzschild solution with p = p = 0, combining the first two 
equations of (fl"4")l implies that 

v' = -A', 

v = -A + z/ (f). (25) 

In this case we absorb the factor i>o(t) into a redinition of the time scale used 
in the metric. Thus we have v = — A and the last two equations of (|14[) become 



2v 

p = —e 



2u' + 1 



p = - ( v" + 2v n + — \ e Zv . (26) 



r 



Note that these two equations both determine the metric function v(r) and the 
pressure (and density) so that one does not a freedom of choice for the pressure. 

We parametrize the exponential metric function v by introducing a potential- 
like function V just as in Eq. (|18l) and in that case the last of the two crucial 
Einstein equations in (1261) become 



p=-(i/' + 2v 12 + ^-) e ^=p= -\(V" + —). (27) 
\ r J 2 r 

Substituting this in the first of Eqs. (|2"6"]l we obtain 

I, .. 2V\ , , /2(1/2)V 1 \ 1 

V V 

= (28) 
This leads to the linear equation 

2V 

V" = (29) 
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Note the difference between this equidimensionai equation and (fit?)) for the ho- 
mogeneous case. This one has the general solution of a linear combination of a 
harmonic oscillator with a Newtonian potential, 

V(r) = kr 2 -—. (30) 
r 

Use of the Schwarzschild radius in this case would impy that we are at distances 
at which the harmonic term is completely negligible. Our metric is thus 

" = -* , — l + J-r./r - 1 "" 

911 = e 2X =e- 2v = l + kr 2 -r s /r = l/g n , 

922 = r 2 = l/g 22 , 

ff33 = r 2 sin 2 6 = l/g 33 . 

g^v = o, \i ^ v. (31) 

This corresponds, without the Newtonian term, to the solution obtained by de 
Sitter for the Einstein equation with a cosmological constant [3] -[5J. 

The pressure and density terms are thus required to be by the equations 
respectively 

1 IV 

V = ~(V" + —) = -3k, 

2 r 

p = -p = 3k. (32) 

Note that the contributions to the energy density and pressure from the New- 
tonian part of the potential vanishes. 

It is of interest that in the nonrelativistic limit this vanishing of the source 
term for that portion of the metric is contrary to what occurs in the Poisson 
equation for a point mass density 

V 2 $ = 4:TTGM5 3 (r), 
GM 

$ = . (33) 

r 

Is there a point mass at the origin in the case of the Einstein equation? This 
would seem to be implied by the Newtonian-Poisson connection. One may be 
tempted to replace V" + 2 ^ r with V 2 V and proclaim that p = —p = —\V 2 V = 
27rr s (5 3 (r)-|-3fc but this is not consistent with the other expression for the pressure 
of p = —p = — = 3k. This calls for a more careful treatment of the Einstein 
equation near the origin. 

2.2 The Schwarzschild Solution as a Valid Approximation 
for a Nonlinear Solution of the Full Non-homogeneous 
Einstein Equation 

In order to treat the problem at the origin more carefully and achieve the aim 
of this paper, we view the Einstein equations in their reduced forms given in 



7 



as the limit for small e of the following 



2v 

p = — e 



2v' 1 \ 1 



p = - [ v" + 2u' 2 + e 2v . (34) 



r 



where 

f= (r 2 +e 2 ) 1/2 . (35) 



We shall solve these equations instead of the reduced forms ([26]) of the Einstein 
equations and view the proper solutions of the Einstein equation as the limit 
of small e of the modified equations. Thus, as before, using Eq. (??) the two 
crucial Einstein equations in (f26f become 



1,„„ . 2V\ 2 J V 1 \ . 1 



-e 



2 V f ' \(1 + V)f f 2 ' '" 2 

V V 



(36) 



and this leads to 



V - ^ _ (37) 
<- (r 2 + e 2 ) 



Clearly, one solution is 
Using the connection 



V 1 {r,s)=k 1 {r 2 +e 2 ). (38) 



V 2 = V 1 I T75) (39) 



dr 
V? 

between the first and second solution of a homogeneous second order differential 
equation, the second solution is 

r dr' 

V 2 {r^)^k 2 {r 2 + e 2 ) J (r , 2+£2)2 . (40) 
We choose the lower limit to be r — oo so that (see Appendix) 



(r + e ) r ^ g / f_^_ 

2e 3 (l + (e/rV 



U 2 (r,e) = ~fc 2 v o Y ; [arctan £ /r- - / ; ]. (41) 



If we let e — > 0, we should get (to match with the Newtonian solution for 
large r) 

, x , a Z" 00 dr' k 2 r 2 k 2 r s 
V 2 (r,0) = -k 2 r 2 =- 2 = -^ = -JL f 42 

J r r 4 3^ 3r r 

so we take 

k 2 = 3r s . (43) 
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(44) 



As a check, let us check that our integrated result (|4"Tj) before we take e — > 
has this same limit 

(r 2 + £ 2 ) e/r 

—3r s — [arctane/r =-) 

2£ 3 (l+(t/r)V 

- ■*.(^ + i)(i-i(i) , -ia-(./D") 

r 2 1 e e 3 e e 3 or 1 " 2 ^ v^ e % 
~ 3r ^ + 2iH" - 3^ - ; + ^3) = -3r -^2i3 + 2^ )( 3^ j 

r 

Thus, our general solution to Eq. (l37t is 

V(r,e) = fci(r 2 + £ 2 ) + 3 ^ (r 2 + £ 2 ) er _ arctan£/r] _ (45) 

with the corresponding metric given by 

911 = e 2X = e- 2 » = l + V(r,e) = l/g n , 

922 = r 2 = l/g 22 , 

333 = r 2 sin 2 # = l/ 5 33 . 

9V = 0, /i/v. (46) 

Now, let us determine what the density and pressure are for the limit of 
small £. This will provide us with insight into the nature of the source terms. 
The simplest way is to evaluate p — — X — ^ . Details omitted here are given 
in the Appendix. There we obtain 

p(r,e) 
= Pi(r,e) +p 2 (r,e) 

= - fc i( (r .2 + g 2) i/2 + !) - 3 ^[ (r . 2+£2) i/2 ^(^ " arctan£/r) + — (p^y - arctan efflfr 

where p\ is the pressure term that arises from the oscillator-like part of the 
solution while p 2 is the pressure term that arises from the Newtonian-like part 
of the solution. For r^O and e — > 

P 2 - -(l- 2 /2r 2 )^+3r s ^3 

-> 0, (48) 
while for r = and £ — ► this becomes 

1 fa} , v 1 r 7T-, 

» - -^-oj-^lo-j] 
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These two limits taken together have the appearance an attractive delta func- 
tion. To complete the verification let us check by integrating the Newtonian 
term for the density, using the divergence theorem, that we obtain the appro- 
priate constant (independent of e) 



J d 3 rp 2 {r,e) = - J tf 



r p 2 (r,e) 



:5, '» / ^[ (,.2^2)1/2 - arctane/r) + —{ -^—^ - arctane/r] 

^[-y-^] = i/^" + -) 

r r z 2 J r 

- [ d 3 r{V" + —) = ! [ d 3 rV 2 V = - lim RHirV'(R) 
2 J r 2 J 2-R-i-co 

= 2tt lim R 2 — ^ arctane/i? = 2tt lim ^ h - — 

R^oo e 3 y R 1 ' R^oo e 3 R R 3 \RJ ' 

= 2irr s = AttGM. (50) 

Such would not be the case for the homogeneous equation which would give 
zero for the integrated density (see below Eq. (|32p). Thus, in the limit e — > 
where our e— modified Einstein equations become the actual Einstein equation, 
the integral of the density over an arbitrarily small volume remains a constant 
independents of e. So, defining (rearranged) 

j3 ( ) = = _ 3 £ / r L^^)_ arctan£/r( L^_ + I)] 7 

v ' 7 4iGM 2^ e 3 L (i + ( e / r ) 2 )i/2 v 2(i + ( £ / r ) 2 )i/2 y 7 V (l + (e/r) 2 )V2 2 7J ' 

(51) 

we have 

r d 3 rS 3 (r,e) = 1. (52) 



Our (5 3 (r,e) therefore has the requisite properties for a distribution that in 
the limit represents a Dirac delta function. Its value for r ^ tends to zero 
as e — > and its integral over all space is unity. This establishes that for 
k± = 0, the source term is non-zero and has the property of a sharply confined 
distribution. What makes this source distinct from others [7], beginning with 
the matching solution of SchwarzschildfS] to an incompressible fluid confined 
within a finite spherical surface, is that within this sharply confined region, the 
pressure is negative. 



2.2.1 The Schwarzschild Limit 

We consider in this section the potential- like function V(r, e) given in Eq. (|45|) . 
We wish to determine in what sense that, if we choose k\ = 0, the second 
portion of V(r, e) for e > sufficiently small agrees with the Schwarzschild 
solution V s (r) = —r s /r for a given range of r. Let us make this statement 
precise. We show that for a positive 5 > that if, 

— =e<6. (53) 
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then 



\V(r,e)-V s (r)\ < 



(54) 



for r in the range e < r < oo. If 5 < 1, then the range for r of agreement between 
the potentials in the above sense would extend down below the Schwarzschild 
radius with no upper bound. 

To show this we consider the Taylor series for V(r, e) — V s (r) in e about 
e = 0. 

V{r,e) - (-)— = 7T^[er- (r 2 + e 2 ) arctane/r] + — (55) 
The series for arctane/r converges for r > e. Expanding we find that 



l/(r )£ )-(-)- 
r 



2 £ 3 [ 1 r 2 ^2n+l iJ J r 

n—0 



3^ 

2 £ 3' 



n=0 n=0 



(-)" /e\ 2 "+ 3 , 



+ 



2e 3L V ^2n + llr/ Z^ 2 n-lVr/ ^ J 



3r s 



n=2 
oo 



n=2 



'(;) 2 g' (2, i + 5, 1 (2„ + 3 ) "- r (?) 2 " 1 

n—0 



Thus, with — = e < <5, we have since e— < 1 



■V(r,e) + 1 



< 



3e 2 



[(£( 



rl L ^(2n + 5)(2n + 3) 



><->"(;) 



2n 



n=0 
2 / oo 



V 1 + D 



(2n + 5)(2n + 3) 
15 



^(2n + 5)(2n + 3) 



/ r* \ 2n 

)(-r(v) ] 
)(-)" (v) 



2 n 



e 2 5 2 

5 <y 



So, in the above sense, we have demonstrated how the Schwarzschild solution 
can be viewed as a valid approximation for a bona fide nonlinear solution of the 
full non-homogeneous Einstein equation. 



3 Discussion 

We have found that the general solution of the Einstein equation for the special 
case of p = —p to yield a metric governed by a linear combination of a New- 
tonian and simple harmonic oscillator potential, two independent solution of a 
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linear second order differential equation for the potential-like function V. If the 
behavior about the origin is not handled carefully, the simple harmonic contri- 
bution gives a constant density and pressure while the Newtonian contribution 
gives rise to no point-like (delta function) source. With k\ — 0, this is the usual 
no source or vacuum solution. Handled more carefully by the method indicated 
we obtain a density and pressure terms sharply peaked about the origin, with 
a unit volume integral. The density corresponding to the elastic portion of the 
metric ranges from k\ for to 3fei. 

Do our mathematical solutions of the inhomogeneous Einstein equation have 
any physical significance? Superimposed on our point-like density is a density 
ranging between two constants, the source of the harmonic oscillator potential- 
like function Vi(r, e) which behaves like k\r 2 for sufficiently small r. The static 
and spherically symmetric metric we started with is distinct from the standard 
time-dependent Friedmann-Lematre-Robertson- Walker (FLRW) metric. It is 
not intended to relate to the universe as a whole but rather to the field produced 
by a single source. There is some superficial similarity between our solution 
and the so-called dark energy and inflation solutions of modern cosmology since 
both involve a pressure with an opposite sign of the density . It may be of 
just academic interest that such potential-like functions not only proportional 
respectively to 1/r but also to the otherwise ubiquitous r 2 potential correspond 
to the solution of the Einstein equation under these circumstances. A positive or 
negative sign of k\ would give a negative or positive pressure and an attractive or 
repulsive force that would increase in magnitude with distance. Obviously since 
there is no evidence for such long range and static increasing forces, k\ must 
be virtually infinitesimal if not zero. Note that in nonrelativistic potential 
theory a constant density would give rise to an attractive Hooke's law force. 
However, that follows from Gauss' law applied to an inverse square field. The 
r 2 potential like function discussed here is completely independent from the 
1/r potential. Another factor to point out is that the functional form of the 
density or pressure due to the harmonic oscillator potential is fixed by our 
solution to the Einstein equations themselves, it is not imposed. The only 
imposition we made on the density and pressure is that they be the negative of 
one another. From a mathematical point of view there is no distinction between 
the solutions discussed here for the inhomogeneous Einstein equation and the 
one we would have obtained by adding a term —pg^ v to the left hand side of the 
Einstein equation and viewing it as an addition to the equation, in analogy to 
the alternative explanation of dark energy. The difference here is that p being 
a constant would be an outcome of the modified Einstein equations and not an 
imposed functional form. 

It was one of Einstein's early goals, although he never succeeded, to incor- 
porate Mach's principle in his general theory of relativity. It has been generally 
regarded that general relativity does not embody Mach's principle, that is that 
geometry can exist independent of matter. It was the Schwarzschild solution 
that seemed to bring this idea its early but reluctant acceptance. That is, a ge- 
ometry arises from the absence of a source term, from the vacuum. Of course, in 
the practical applications of the Schwarzschild solution to the precession prob- 
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lem of Mercury and the bending of light, it was always assumed that looming 
behind the formal sourceless equation was a real sun. Nevertheless, a possible 
formal interpretation is that a curved space exists without an identifiable source, 
thus obviating the need for Mach's principle 

Our result has been to replace the Schwarzschild solution to the sourceless 
spherically symmetric static environment, which then, as now seems to allow 
the existence of non-trivial spacetime curvature in absence of any matter, with 
a solution that does not correspond to a sourceless environment but yet leads 
nevertheless to a metric that can approach the Schwarzschild with arbitrary 
accuracy in an asymptotic way. In doing so, for this particular case at least, 
Mach's principle, the idea that geometry emerges as an interaction between an 
identifiable matter term and geometry is preserved[9]. 

A Details of Solution V^r, e) of e— Modified Ein- 
stein Equations 

Here we consider some details omitted from the text. The expression for V2 (r, e) 
is obtained below, 



V 2 (r,e) = = 
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Next we combine terms to determine the expression for the pressure given in 
Eq. (??) 



P 



,-2 



V V 



V , , 1 . er 

k\ + k 2 — 7-5- — - arctane/r , 
2e 3 (r 2 + £ ) 

1 „, , 2r r er . , , (r 2 + e 2 ) , e 2er 2 

V' = 2fcir + fc 2 — [ -arctane/r] +fc 2 - 



'2e 3L (r 2 +e 2 ) ' J 2e 3 L (r 2 +e 2 ) (r 2 + e 2 ) 2 (r 2 + e 2 ) J 

r er (r 2 + e 2 ) e 2 

2fcir + k2-r[-r-n 2T ~ arctane/r] + k 2 1 1 



e 3 (r 2 +e 2 ) e 2 L (r 2 +e 2 ) 2 

„, r , er , . k 2 

2kir + k 2 ^(-r^ - arctane/r) + 



; e 3V (r 2 +e 2 ) ' ; (r 2 +e 2 ) 

r 2 /e 2 r 1 

7-5 ^7 r arctane/r + — -= 

' (r 2 + e 2 ) e 3 (r 2 + e 2 ) 



= 2kir + k 2 (— arctane/r) = 2fcir H —( arctane/r) (A. 2) 

e e 3 e 3 r 
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